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Abstract 

We present applications of variational-wavelet approach to nonlinear (ra- 
tional) rms envelope equations. We have the solution as a multiresolution 
(multiscales) expansion in the base of compactly supported wavelet basis. We 
give extension of our results to the cases of periodic beam motion and arbi- 
trary variable coefficients. Also we consider more flexible variational method 
which is based on biorthogonal wavelet approach. 
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1 Introduction 



In this paper we consider the applications of a new numerical-analytical tech- 
nique which is based on the methods of local nonlinear Fourier analysis or 
Wavelet analysis to the nonlinear beam/accelerator physics problems related 
to root-mean-square (rms) envelope dynamics [1]. Such approach may be use- 
ful in all models in which it is possible and reasonable to reduce all complicated 
problems related with statistical distributions to the problems described by 
systems of nonlinear ordinary /partial differential equations. In this paper we 
consider approach based on the second moments of the distribution functions 
for the calculation of evolution of rms envelope of a beam. 

The rms envelope equations are the most useful for analysis of the beam 
self-forces (space-charge) effects and also allow to consider both transverse 
and longitudinal dynamics of space-charge-dominated relativistic high-bright 
ness axisymmetric/asymmetric beams, which under short laser pulse-driven 
radio-frequency photoinjectors have fast transition from nonrelativistic to rel- 
ativistic regime [2]- [3]. 

From the formal point of view we may consider rms envelope equations 
after straightforward transformations to standard Cauchy form as a system of 
nonlinear differential equations which are not more than rational (in dynam- 
ical variables). Such rational type of nonlinearities allow us to consider some 
extension of results from [4]-[12], which are based on application of wavelet 
analysis technique to variational formulation of initial nonlinear problem. 

Wavelet analysis is a relatively novel set of mathematical methods, which 
gives us a possibility to work with well-localized bases in functional spaces and 
give for the general type of operators (differential, integral, pseudodifferential) 
in such bases the maximum sparse forms. 

An example of such type of basis is demonstrated on Fig. 1. 

Our approach in this paper is based on the generalization [13] of variatio- 
nal- wavelet approach from [4]- [12], which allows us to consider not only poly- 
nomial but rational type of nonlinearities. 

So, our variational-multiresolution approach gives us possibility to con- 
struct explicit numerical-analytical solution for the following systems of non- 
linear differential equations 

i = R(z, t) or Q(z, t)z = P{z, t), (1) 

where z(t) = (zi(t), ...,z n (t)) is the vector of dynamical variables Zi(t), 
R(z, t) is not more than rational function of z, 

P(z,t),Q(z,t) are not more than polynomial functions of z and P,Q,R 
have arbitrary dependence of time. 
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The solution has the following form 

z(t)=z$r(t) + j £*iM, "i-y ( 2 ) 

j>N 

which corresponds to the full multiresolution expansion in all time scales. 
Formula (g) gives us expansion into a slow part z s Jf w and fast oscillating 
parts for arbitrary N. So, we may move from coarse scales of resolution to 
the finest one for obtaining more detailed information about our dynamical 
process. The first term in the RHS of equation (2) corresponds on the global 
level of function space decomposition to resolution space and the second one 
to detail space. In this way we give contribution to our full solution from each 
scale of resolution or each time scale (detailed description we give in part 
3.2 and numerical illustration in part 7 below). The same is correct for the 
contribution to power spectral density (energy spectrum): we can take into 
account contributions from each level/scale of resolution. 

In part 2 we describe the different forms of rms equations. Starting in part 
3.1 from variational formulation of initial dynamical problem we construct via 
multiresolution analysis (3.2) explicit representation for all dynamical vari- 
ables in the base of compactly supported (Daubechies) wavelets. Our solu- 
tions (3.3) are parametrized by solutions of a number of reduced algebraical 
problems one from which is nonlinear with the same degree of nonlinearity 
and the rest are the linear problems which correspond to particular method 
of calculation of scalar products of functions from wavelet bases and their 
derivatives. Then we consider further extension of our previous results. In 
part 4 we consider modification of our construction to the periodic case, in 
part 5 we consider generalization of our approach to variational formulation 
in the biorthogonal bases of compactly supported wavelets and in part 6 to 
the case of variable coefficients. In part 7 we consider results of numerical 
calculations. 

2 RMS Equations 

Below we consider a number of different forms of RMS envelope equations, 
which are from the formal point of view not more than nonlinear differential 
equations with rational nonlinearities and variable coefficients. Let f(x\,X2) 
be the distribution function which gives full information about noninteract- 
ing ensemble of beam particles regarding to trace space or transverse phase 
coordinates (x\,X2). Then (n,m) moments are: 

xlxff(x 1 ,x 2 )&x 1 &X2 (3) 
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Figure 1: Wavelets at different scales and locations. 



The (0,0) moment gives normalization condition on the distribution. The 
(1,0) and (0,1) moments vanish when a beam is aligned to its axis. Then we 
may extract the first nontrivial bit of 'dynamical information' from the second 
moments 



a x 1 = 

2 



<x?>=//x?/(x 1 ,x 2 )d, 1 dx 2 
x l >= J J xlf(xi,x 2 )dx 1 dx 2 
xix 2 >= J j Xix 2 f(x\,x 2 )dxidx 2 
RMS emittance ellipse is given by 



X1X2 



< 
< 



(4) 



< xl X xl > 



< X\X 2 >' 



(5) 



Expressions for twiss parameters are also based on the second moments. 

We will consider the following particular cases of rms envelope equations, 
which described evolution of the moments (4) ([l]-[3] for full designation): 

for asymmetric beams we have the system of two envelope equations of 



the second order for a Xl and a X2 : 
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the envelope equation for an axisymmetric beam is 

" . 'i , n 2 { "f'V k s , £ n,th 

cr + a 1- S2 — cr = ^ + -^-^ (7) 

7 V 7 / °"7 <t j 7^ 

Also we have related Lawson's equation for evolution of the rms envelope in 
the paraxial limit, which governs evolution of cylindrical symmetric envelope 
under external linear focusing channel of strenghts K r : 

where 

K r = -F r /r(3 2 -fmc 2 , (3 = u b /c = ^l-l~ 2 (9) 

After transformations to Cauchy form we can see that all this equations 
from the formal point of view are not more than ordinary differential equa- 
tions with rational nonlinearities and variable coefficients and correspond to 
the form (1) (also, we may consider regimes in which 7, 7' are not fixed func- 
tions/constants but satisfy some additional differential constraint/equation, 
but this case does not change our general approach). 

3 Rational Dynamics 

The first main part of our consideration is some variational approach to this 
problem, which reduces initial problem to the problem of solution of functional 
equations at the first stage and some algebraical problems at the second stage. 
We have the solution in a compactly supported wavelet basis. Multiresolu- 
tion expansion is the second main part of our construction. The solution is 
parameterized by solutions of two reduced algebraical problems, one is nonlin- 
ear and the second are some linear problems, which are obtained from one of 
the next wavelet constructions: the method of Connection Coefficients (CC), 
Stationary Subdivision Schemes (SSS). 

3.1 Variational Method 

Our problems may be formulated as the systems of ordinary differential equa- 
tions 

Qi(x)^ = Pi(x,t), x = (xi,...,x n ), (10) 
i = l,...,n, maxdeg Pi = p, rn.ax.deg Qi = q 

i i 
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with fixed initial conditions Xj(0), where Pi, Qi are not more than polynomial 
functions of dynamical variables Xj and have arbitrary dependence of time. 
Because of time dilation we can consider only next time interval: < t < 1. 
Let us consider a set of functions 

= Xi^iQiVi) + Piyi (11) 

and a set of functionals 

Fi{x) = I 1 <$>i{t)dt - QiXiyi \l, (12) 
Jo 

where Ui(t) (yi(0) = 0) are dual (variational) variables. It is obvious that the 
initial system and the system 

Fiix) = (13) 

are equivalent. Of course, we consider such Qi(x) which do not lead to the 
singular problem with Qi(x), when t = or t = 1, i.e. Qi(x(0)), Qi(x(l)) ^ oo. 

In part 5 we consider more general approach, which is based on possibility 
taking into account underlying symplectic structure and on more useful and 
flexible analytical approach, related to bilinear structure of initial functional. 
Now we consider formal expansions for Xi,yf. 

Xiit) = Xi (Q) Vj{t) = ^2vjMt), (14) 

k r 

where ip^ (t) are useful basis functions of some functional space {I? , D\ Sobolev, 
etc) corresponding to concrete problem and because of initial conditions we 
need only (/Jfc(0) = 0. 

A = {A,} = {Aa = (A, 1 ,A 4 2 ,...,Af), r = l,...,N, i = l,...,n, (15) 

where the lower index i corresponds to expansion of dynamical variable with 
index i, i.e. X{ and the upper index r corresponds to the numbers of terms 
in the expansion of dynamical variables in the formal series. Then we put 
(14) into the functional equations (jl^) and as result we have the following 
reduced algebraical system of equations on the set of unknown coefficients 
of expansions (|l4 



L(Q ij ,X,a I ) = M(P ij ,X,p J ), (16) 
where operators L and M are algebraization of RHS and LHS of initial problem 



(|T0|), where A (15) are unknowns of reduced system of algebraical equations 
(RSAE)(||). 
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Qij are coefficients (with possible time dependence) of LHS of initial sys- 
tem of differential equations ( |l0|) and as consequence are coefficients of RSAE. 

Pij are coefficients (with possible time dependence) of RHS of initial sys- 
tem of differential equations ( |l0|) and as consequence are coefficients of RSAE. 

/ = (ii, ...,i q+2 ), J = (ji, jp+i) are multiindexes, by which are labelled 
a j and /3j — other coefficients of RSAE ([u]): 

Pj = {Pn... jp+ A = f II <Ph> ( 17 ) 
l<jk<p+l 

where p is the degree of polinomial operator P (fil 



ai = {a h ...a iq+2 } = J Vh-Vis-Pi^, ( 18 ) 

h,— Aq+2 

where q is the degree of polynomial operator Q (|T(i|), %i = (1, ...,q + 2), <p\ s = 
difijdt. 

Now, when we solve RSAE (0) and det ermine unknown coefficients from 
formal expansion ( |l~4|) we therefore obtain the solution of our initial problem. 
It should be noted if we consider only truncated expansion (|l4]) with N terms 
then we have from (]l~6| ) the system of JVx n algebraical equations with degree 
£ = max{p, q} and the degree of this algebraical system coincides with degree 
of initial differential system. So, we have the solution of the initial nonlinear 
(rational) problem in the form 

N 

x i (t)=x i (Q) + y t ^X k (t), (19) 



fc=i 



where coefficients are roots of the corresponding reduced algebraical (poly- 
nomial) problem RSAE (|l^). Consequently, we have a parametrization of so- 
lution of initial problem by solution of reduced algebraical problem ( |l6| ) . The 
first main problem is a problem of computations of coefficients a/ (|l8|), 0j 
(17) of reduced algebraical system. As we will see, these problems may be 



explicitly solved in wavelet approach. 

Next we consider the construction of explicit time solution for our problem. 
The obtained solutions are given in the form (|l9|), where X^(t) are basis func- 
tions and \\ are roots of reduced system of equations. In our first wavelet 
case Xk(t) are obtained via multiresolution expansions and represented by 
compactly supported wavelets and A^. are the roots of corresponding general 
polynomial system ( |i~6| ) with coefficients, which are given by CC or SSS con- 
structions. According to the variational method to give the reduction from 
differential to algebraical system of equations we need compute the objects 
a / and (3j. 
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3.2 Wavelet Framework 



Our constructions are based on multiresolution approach. Because affine 
group of translation and dilations is inside the approach, this method re- 
sembles the action of a microscope. We have contribution to final result from 
each scale of resolution from the whole infinite scale of spaces. More exactly, 
the closed subspace Vj(j G Z) corresponds to level j of resolution, or to scale 
j. We consider a r-regular multiresolution analysis of L 2 (R n ) (of course, we 
may consider any different functional space) which is a sequence of increasing 
closed subspaces Vji 

...V-2 C V-i C Vq C Vl C V 2 C ... (20) 
satisfying the following properties: 

jez jez 

f{x) G Vj <=> f(2x) G V j+1 , 

f{x) G V < = > f(x -k)€ V , ,Vk G Z" (21) 

There exists a function ip G Vq such that {<^ ,fc(x) = (fix — k),k G Z™} forms 
a Riesz basis for Vq. 

The function ip is regular and localized: ip is C r ~ 1 , ip( r ~^ is almost 
everywhere differentiable and for almost every x G R n , for every integer a < r 
and for all integer p there exists constant C p such that 

I d a pix) \< C p (l + \x\)- p (22) 

Let (pix) be a scaling function, ipix) is a wavelet function and <Piix) = 
pix — i). Scaling relations that define ip, ip are 

N-l N-l 

= a wi^ x - k ) = a k^ki2x), (23) 

N-2 

1>ix) = i-l) k a k +M2x + k). (24) 

k=-i 

Let indices £, j represent translation and scaling, respectively and 

(Pjlix) =2 j / 2 <pi2 j x-£) (25) 

then the set {pj t k}, k G Z n forms a Riesz basis for Vj. The wavelet function ip 
is used to encode the details between two successive levels of approximation. 
Let Wj be the orthonormal complement of Vj with respect to V^+i: 

V j+ i = Vj®Wj. (26) 
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Then just as Vj is spanned by dilation and translations of the scaling function, 
so are Wj spanned by translations and dilation of the mother wavelet ipjk(x), 
where 

i>jk{x) = 2 j / 2 ip(2 j x - k). (27) 
All expansions which we used are based on the following properties: 

{VVfc}; 3ik £ Z is a Hilbertian basis of L 2 (R) 
{<Pjk}j>o,kez is an orthonormal basis forL 2 (R), 

oo 

L 2 (R) = F O 0^, (28) 
or 

Wo,ki V'j,fe}j>o,feez is an orthonormal basis forL 2 (R). 

3.3 Wavelet Computations 

Now we give construction for computations of objects (17), (18) in the wavelet 
case. We use compactly supported wavelet basis: orthonormal basis for func- 
tions in L 2 (R). 

Let be / : R — ► C and the wavelet expansion is 

oo 

f( x ) = + c jktpjk{x) (29) 

fez j=okez 

If in formulae (29) Cjk = for j > J, then f(x) has an alternative ex- 
pansion in terms of dilated scaling functions only f{x)= J2 cj£(pji(x). This 

is a finite wavelet expansion, it can be written solely in terms of translated 
scaling functions. Also we have the shortest possible support: scaling function 
DN (where N is even integer) will have support [0, N — 1] and N/2 vanishing 
moments. There exists A > such that DN has XN continuous derivatives; 
for small N, A > 0.55. To solve our second associated linear problem we need 
to evaluate derivatives of f(x) in terms of ip(x). Let be ip™ = d n (pe(x)/dx n . 
We consider computation of the wavelet - Galerkin integrals. Let f d (x) be 
d-derivative of function f(x), then we have f d (x) = 2~2e c lPe( x )i anc ^ va l ues 
<p d {x) can be expanded in terms of <p(x) 

4{x) = J2 X ^ m (x), (30) 

m 
oo 

A m = J <p d (x)<p m (x)dx, 

— oo 
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where \ m are wavelet-Galerkin integrals. The coefficients X m are 2-term con- 
nection coefficients. In general we need to find (di > 0) 

oo 

Ktt= J IW;(^ ( 31 ) 

— oo 

For Riccati case we need to evaluate two and three connection coefficients 

oo 

/oo r 
<p dl (x)<pf(x)dx, A dld2d3 = <f dl (x)(pf 2 {x)^(x)dx (32) 
-oo J 

— oo 

According to CC method [14] we use the next construction. When ./V in scaling 
equation is a finite even positive integer the function (p(x) has compact support 
contained in [0,N — 1]. For a fixed triple (di,d 2 ,ds) only some A^ 2 * 1 are 
nonzero: 2 - N < I < N - 2, 2-N<m<N-2, \t - m\ < N - 2. There 
are M = 3N 2 - 9N + 7 such pairs (£, m). Let A dld2d ' 3 be an M- vector, whose 
components are numbers A^ 2 ^ 3 . Then we have the first reduced algebraical 
system : A satisfy the system of equations (d = d\ + d2 + d%) 

AK d ld2 d, = 2 i-d A d ld2 d 3 f = J2 a p a q , 2e+p a r ^ m+p (33) 

v 

By moment equations we have created a system of M + d + 1 equations in M 
unknowns. It has rank M and we can obtain unique solution by combination 
of LU decomposition and QR algorithm. The second reduced algebraical 
system gives us the 2-term connection coefficients: 

Ah d ld2 = 2 i-d A dKh ) d = dl + d2: A e , q = Y,a p a q - 2 e +p (34) 

p 

For nonquadratic case we have analogously additional linear problems for ob- 
jects (31). Solving these linear problems we obtain the coefficients of nonlinear 
algebraical system (16) and after that we obtain the coefficients of wavelet ex- 
pansion (19). As a result we obtained the explicit time solution of our problem 
in the base of compactly supported wavelets. We use for modelling D6, D8, 
D10 functions and programs RADAU and DOPRI for testing. 

In the following we consider extension of this approach to the case of peri- 
odic boundary conditions, the case of presence of arbitrary variable coefficients 
and more flexible biorthogonal wavelet approach. 
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4 Variational Wavelet Approach for Peri- 
odic Trajectories 

We start with extension of our approach to the case of periodic trajectories. 
The equations of motion corresponding to our problems may be formulated 
as a particular case of the general system of ordinary differential equations 
dxi/dt = fi(xj,t), = l,...,n), < t < 1, where fi are not more than 
rational functions of dynamical variables Xj and have arbitrary dependence 
of time but with periodic boundary conditions. According to our variational 
approach we have the solution in the following form 

Xi(t) = Xi(0) + J2 X iMt), Xi(0) = Xi(l), (35) 
k 

where are again the roots of reduced algebraical systems of equations with 
the same degree of nonlinearity and ipk(i) corresponds to useful type of wavelet 
bases (frames). It should be noted that coefficients of reduced algebraical 
system are the solutions of additional linear problem and also depend on 
particular type of wavelet construction and type of bases. 

This linear problem is our second reduced algebraical problem. We need 
to find in general situation objects 

oo 

Ktt= / IW;(*) d *> ( 36 ) 

— oo 

but now in the case of periodic boundary conditions. Now we consider the 
procedure of their calculations in case of periodic boundary conditions in 
the base of periodic wavelet functions on the interval [0,1] and corresponding 
expansion (35) inside our variational approach. Periodization procedure gives 
us 

<Pj,k( x ) = ^ZwiA 00 ~ ( 37 ) 

eez 

So, <p,ip are periodic functions on the interval [0,1]. Because <pj t k = <fj,k' if 
k = fc'mod^), we may consider only < k < 2 3 and as consequence our 
multiresolution has the form Vj = L 2 [0, 1] with Vj = sp&n{tf>j^} 2 k=o [15]. 

i>o 
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Integration by parts and periodicity gives useful relations between objects 
(36) in particular quadratic case (d = d± + d-i): 

A fci,fc 2 - y~ L ) A k u k 2 i A k u k 2 - A 0,fe 2 -fci = A k 2 -k ± W»J 

So, any 2-tuple can be represent by A^. Then our second additional linear 
problem is reduced to the eigenvalue problem for {A^} <fc< 2 j by creating a 
system of 2- 7 homogeneous relations in A^ and inhomogeneous equations. So, 
if we have dilation equation in the form ip(x) = \[2 J2kez hk<fi{2x — k), then 
we have the following homogeneous relations 

N-l N-l 

*4 = 2 " E E h m h e A d £+2k _ m , (39) 

m=0 £=Q 

or in such form AX d = 2 d X d , where X d = {A. d } 0<k<2 j ■ Inhomogeneous equa- 
tions are: 

Y / MfAf = d\2- j / 2 , (40) 

t 

where objects Mf(\i\ < N - 2) can be computed by recursive procedure 

Mf = 2-^ 2d+1 )/ 2 M/, Ml =< x k , <p 0t e >= £ ( k ) n k ~ j Ml = 1. 

j=o V/ 

(41) 

So, we reduced our last problem to standard linear algebraical problem. Then 
we use the same methods as in part 3.3. As a result we obtained for closed 
trajectories of orbital dynamics the explicit time solution (35) in the base of 
periodized wavelets (37). 



5 Variational Approach in Biorthogonal 
Wavelet Bases 

Now we consider further generalization of our variational wavelet approach. 
Because integrand of variational functionals is represented by bilinear form 
(scalar product) it seems more reasonable to consider wavelet constructions 
[16] which take into account all advantages of this structure. 
The action functional for loops in the phase space is 

F( 7 )= / ' pdq- f 1 H(t,j(t))dt (42) 
J7 Jo 
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The critical points of F are those loops 7, which solve the Hamiltonian equa- 
tions associated with the Hamiltonian H and hence are periodic orbits. 

Let us consider the loop space = C°°(S 1 , R 2n ), where S 1 = R/Z, of 
smooth loops in R 2n . Let us define a function $ : Q — > R by setting 

(f>(x) = f 1 - < -Jx,x>dt- f 1 H(x(t))dt, (43) 
Jo 2 Jo 

The critical points of $ are the periodic solutions of x = Xh(x). Computing 
the derivative at x G £1 in the direction of y £ fi, we find 

d r 1 

&(x)(v) = —Hx + ey)\ e=0 = / < -Jx-vH{x),y> dt (44) 
oe jo 

Consequently, $'(x)(y) = for all y € iff the loop x satisfies the equation 

- Jx(t) - s?H(x(t)) = 0, (45) 

i.e. x(t) is a solution of the Hamiltonian equations, which also satisfies x(0) = 
x(l), i.e. periodic of period 1. 

But now we need to take into account underlying bilinear structure via 
wavelets. 

We started with two hierarchical sequences of approximations spaces [16]: 

. . . V-2 C V-x C V C Vi C V 2 . . . , . . . V-2 C V-! C V C V x C V 2 ■ ■ ■ , 

and as usually, Wq is complement to Vq in Vi, but now not necessarily or- 
thogonal complement. New orthogonality conditions have now the following 
form: 

W ±Vq, W ±V , Vj±Wj, Vj^Wj (46) 

translates of tp span Wq, translates of ip span Wo. Biorthogonality condi- 
tions are 

/oo 
ipjk(x)4> j 'k'(x)dx = 6 kkl 6 j:jl , (47) 
-00 

where ipjk( x ) = 2^ 2 i/}(2 : >x — k). Functions tp(x),<p(x — k) form dual pair: 

<<p(x-k),0(x-£)>=5id, < (p(x-k),4>{x-£) >= for VA;, VI (48) 

Functions (p,ip generate a multiresolution analysis. (p(x — k), ip(x — k) are 
synthesis functions, (p(x — £), ip(x — £) are analysis functions. Synthesis func- 
tions are biorthogonal to analysis functions. Scaling spaces are orthogonal 
to dual wavelet spaces. Two multiresolutions are intertwining Vj + Wj = 
Vj+i, Vj + Wj = Vj + \. These are direct sums but not orthogonal sums. 
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So, our representation for solution has now the form 

f(t) = J2hk^jk(t), (49) 

where synthesis wavelets are used to synthesize the function. But bjk come 
from inner products with analysis wavelets. Biorthogonality yields 

km = J f(t)j> im (t)dt. (50) 

So, now we can introduce this more complicated construction into our vari- 
ational approach. We have modification only on the level of computing co- 
efficients of reduced nonlinear algebraical system. This new construction is 
more flexible. Biorthogonal point of view is more stable under the action 
of large class of operators while orthogonal (one scale for multiresolution) is 
fragile, all computations are much more simpler and we accelerate the rate 
of convergence. In all types of (Hamiltonian) calculation, which are based on 
some bilinear structures (symplectic or Poissonian structures, bilinear form of 
integrand in variational integral) this framework leads to greater success. 



6 Variable Coefficients 

In the case when we have situation when our problem is described by a system 
of nonlinear (rational) differential equations, we need to consider extension 
of our previous approach which can take into account any type of variable 
coefficients (periodic, regular or singular). We can produce such approach 
if we add in our construction additional refinement equation, which encoded 
all information about variable coefficients [17]. According to our variational 
approach we need to compute integrals of the form 

f b ij (t)(<p 1 ) d i(2 m t-k 1 )(<p 2 ) d *(2 m t-k 2 )dx, (51) 

JD 

where now hj(t) are arbitrary functions of time, where trial functions pi,<p2 
satisfy a refinement equations: 

<Pi(t) = OiM(2t - k) (52) 
fcez 

If we consider all computations in the class of compactly supported wavelets 
then only a finite number of coefficients do not vanish. To approximate the 



14 



non-constant coefficients, we need choose a different refinable function 
along with some local approximation scheme 



(B e f)(x) := £ F etk (f)<p3(2tt - k), (53) 

where & are suitable functionals supported in a small neighborhood of 2~^k 
and then replace by in ( |5l| ) by Bgbij(t). In particular case one can take 
a characteristic function and can thus approximate non-smooth coefficients 
locally. To guarantee sufficient accuracy of the resulting approximation to 



(51) it is important to have the flexibility of choosing tp^ different from tpx,ip2- 



In the case when D is some domain, we can write 

bij(t)\o= £ b lj (t) XD (2 e t-k), (54) 

0<k<2 e 

where xd is characteristic function of D. So, if we take = xd, which is 



again a refinable function, then the problem of computation of (51) is reduced 
to the problem of calculation of integral 

H{k 1 ,k 2 ,k 3 ,k 4 ) = H{k)= (55) 
(p 4 (2 j t - ki)<p 3 (2 e t - k 2 )^ dl {2 r t - k z )^ d 2 2 {2 s t - & 4 )dx 

The key point is that these integrals also satisfy some sort of refinement 
equation [17]: 

2-^H(k) = J2b2k-eH(£), / x = d 1 + cZ 2 . (56) 
eez 

This equation can be interpreted as the problem of computing an eigen- 
vector. Thus, we reduced the problem of extension of our method to the case 
of variable coefficients to the same standard algebraical problem as in the 
preceding sections. So, the general scheme is the same one and we have only 
one more additional linear algebraic problem by which we in the same way 
can parameterize the solutions of corresponding problem. 



7 Numerical Calculations 

In this part we consider numerical illustrations of previous analytical ap- 
proach. Our numerical calculations are based on compactly supported Daubechies 
wavelets and related wavelet families. 
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On Fig. 2 we present according to formulae (2) contributions to approxi- 
mation of our dynamical evolution (top row on the Fig. 3) starting from the 
coarse approximation, corresponding to scale 2° (bottom row) to the finest 
one corresponding to the scales from 2 1 to 2 5 or from slow to fast components 
(5 frequencies) as details for approximation. Then on Fig. 3, from bottom 
to top, we demonstrate the summation of contributions from corresponding 
levels of resolution given on Fig. 2 and as result we restore via 5 scales (fre- 
quencies) approximation our dynamical process (top row on Fig. 3 ). In this 
particular model case we considered for approximation simple two frequencies 
harmonic process. But the same situation we have on the Fig. 5 and Fig. 6 
in case when we added to previous 2-frequencies harmonic process the noise 
as perturbation. Again, our dynamical process under investigation (top row 
of Fig. 6) is recovered via 5 scales contributions (Fig. 5) to approximations 
(Fig. 6). The same decomposition/approximation we produce also on the level 
of power spectral density in the process without noise (Fig. 4) and with noise 
(Fig. 7). On Fig. 8 we demonstrate the family of localized contributions to 
beam motion, which we also may consider for such type of approximation. 

It should be noted that complexity of such algorithms are minimal regard- 
ing other possible. Of course, we may use different multiresolution analysis 
schemes, which are based on different families of generating wavelets and apply 
such schemes of numerical-analytical calculations to any dynamical process 
which may be described by systems of ordinary /partial differential equations 
with rational nonlinearities [13]. 
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Figure 2: Contributions to approximation: from scale 2 1 to 2 5 (without noise). 
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Figure 4: Power spectral density: from scale 2 1 to 2 5 (without noise) 
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Figure 6: Approximations: from scale 2 1 to 2 5 (with noise). 
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Figure 7: Power spectral density: from scale 2 1 to 2 5 (with noise) 
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Figure 8: Localized contributions to beam motion. 
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